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Abstract: We show that for a class of two-loop diagrams, the on-shell part of the integration- 
by-parts (IBP) relations correspond to exact meromorphic one-forms on algebraic curves. 
Since it is easy to find such exact meromorphic one-forms from algebraic geometry, this 
idea provides a new highly efficient algorithm for integral reduction. We demonstrate the 
power of this method via several complicated two-loop diagrams with internal massive legs. 
No explicit elliptic or hyperelliptic integral computation is needed for our method. 
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1 Introduction 

With the beginning of Run II of the Large Hadron Collider (LHC), we need high precision 
scattering amplitudes in Quantum Chromodynamics and the Standard Model, to reduce 
the theoretical uncertainty. The precise scattering amplitude computation suffers from 
problems of the large number of loop Feynman diagrams and difficult integrations for each 
loop diagram. This paper aims at developing a new method of reducing loop integrals to 
the minimal set of integrals, i.e., master integrals (Mis). 

Traditionally, integral reduction can be achieved by applying integration-by-parts (IBP) 
identities [1] and considering the other symmetries of loop diagrams. However, given a two- 
loop or higher-loop integral, it is difficult to find a particular IBP identity which reduce 
it to Mis without introducing unwanted terms. There are several implements of IBPs 
generating codes AIR [2], FIRE [3-5] and Reduze [6, 7], based on Laporta algorithm [8], by 
the computation of Gaussian elimination or Grobner basis. For multi-loop diagrams with 
high multiplicities or many mass scales, it may take a lot of time and computer RAM to 
finish the integral reduction. There are also several new approaches for integral reduction, 
based on the study of the Lie algebra structure of IBPs [9], Syzygy computation [10, 11], 
reductions over finite fields [12], and differential geometry [13]. Besides, the number of 
master integrals can be determined by the critical points of polynomials [14]. 
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We present a new method of integral reduction, for a class of multi-loop diagrams, 
based on unitarity [15-32] and the analysis of algebraic curves [33-35]. We show that for 
a ID-dimensional L-loop diagram, if the unitarity cut solution V is an irreducible algebraic 
curve, then the on-shell IBPs of this diagram correspond to exact meromorphic 1-forms on 
V. For an algebraic curve, it is very easy to find the exact meromorphic 1-froms, based 
on algebraic geometry. Hence for this class of diagrams, we can derive the on-shell part of 
IBPs very efficiently from our method. 

Schematically, an IBP relation. 


f d^h 

J {27r)D- 


d^li d uf 

(27r)^ ^ 


( 1 . 1 ) 


on the unitarity cut V : Di = ... = Dk = 0 becomes contour integral relations [25-32], 


j)u = 0, ( 1 . 2 ) 

where the integrals are along combinations of non-trivial cycles of V, contours surrounding 
singular points of V and poles of lv. If V is an algebraic curve, then the contours are 
one-dimensional. Furthermore, if V is irreducible, then V has a complex structure and tv 
is a meromorphic 1-form [33]. In this case, (1.2) holds for all contours, so it implies the oj 
is an exact form, 

u! = dF, (1.3) 

where F is a meromorphic function on V. So the on-shell part of IBP relations for this 
diagram correspond to exact meromorphic 1-forms. Mathematically, it is very easy to find 
exact meromorphic 1-forms on an algebraic curve, so we can quickly get the on-shell IBPs 
for this diagram. 

In this paper, we show several two-loop examples with internal masses for our new 
method. Multi-loop integrals with internal masses appear frequently in QCD/SM scatter¬ 
ing amplitudes, and are bottlenecks for integral reduction or evaluation. We use these 
complicated cases to show the power of our method: 


1. D = A planar double box with internal massive legs. The unitarity cut for this diagram 
is an elliptic curve. The maximal unitarity structure of the symmetric double box, 
with internal massive legs, was studied in [31]. Here we derive the integral reduction 
for the general cases, based on the analysis of differential forms on an elliptic curve. 
We also reduce integrals with doubled propagators based on algebraic curves, which 
were not considered in [31]. 

2. D = 2 sunset diagram. The unitarity cut for this diagram is again an elliptic curve. 
We derive the analytic integral reduction based on the analysis of elliptic curves. 

3. D = A non-planar crossed box with internal massive legs. The unitarity cut for this 
diagram is a genus-3 hyperelliptic curve. The unitarity cut is more complicated than 
the planar two-loop counterparts, however, our method also works for this case. We 
get the analytic integral reduction from the analysis of the hyperelliptic curve. 
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In these examples, we get all the on-shell IBPs analytically. The algorithm is realized by 
a Mathematica code containing algebraic geometry tools. For each diagram, the analytic 
integral reduction is extremely fast, which has the time order of minutes. 

We have the following remarks, 

• Although the mathematic objects are elliptic or hyperelliptic, we do not need the 
explicit form of elliptic/hyperelliptic functions, or elliptic/hyperelliptic integrals. Only 
the differential relations for elliptic and hyperelliptic functions are needed. These 
relations involve rational coefficients only and are easy to find. 

• The method presented in this paper is different from the maximal unitarity method. 
For the maximal unitarity method, we need to perform contour integrals to extract 
the master integral coefficients. Our method use the integrand reduction [36-39] via 
Grobner basis [39-45] first, to reduce the loop amplitude to an integrand basis. Then 
we use the knowledge of algebraic curves, to reduce the integrand basis further to 
master integrals, fn this way, we avoid the explicit elliptic or hyperelliptic integral 
computations. 

This paper is organized as follows: fn section 2, we present our method based on 
algebraic curves. In section 3 and 4, the double box diagram (elliptic) and sunset diagram 
(elliptic) with internal masses will be explicitly presented. In section 5, we consider the 
integral reduction for the massive nonplanar box diagram (hyperelliptic). The rudiments 
of the knowledge of algebraic curves are included in the appendix. 

2 Integral Reduction via the Analysis of Algebraic Curves 

Generically, for a quantum field theory, the L-loop amplitude can be written as [15, 16], 


loop _ Cklk + rational terms 


( 2 . 1 ) 


k 


The set {Ik} is called the master integral (MI) basis whose elements are independent loop 
integrals. In practice, for amplitudes with multiple loops, high multiplicities or several mass 


scales, it is quite difficult to determine the set of master integral or reduce a generic integral. 


5 

to the linear combination of master integrals. 

Traditionally, the integral reduction is done by using IBP identities [1] 



( 2 . 2 ) 



(2.3) 


if there is no boundary term, fn general, it is difficult to find the IBP relations for a 
multi-loop integral reduction. 






We present a new way of integral reduction, based on maximal unitarity method and 
algebraic curves. Given a Feynman integral with k propagators, maximal unitarity method 
split (2.1) as [17-32] , 

Int = Cili + (integrals with fewer-than-k propagators) + rational terms (2-4) 
i 

where the first sum is over the master integral with exact k propagators. 

The condition the all internal legs are on-shell, is called the maximal unitarity cut, 

F ; Di = ... = Dk = Q, (2.5) 


and the solution set for this equation system is an algebraic variety V. V can be a set 
of discrete points, algebraic curves or surfaces. (See [46, 47] for the detailed mathematical 
study of multi-loop unitarity cut solutions.) Maximal unitarity replaces the original integral 
with contour integrals [25-32], schematically, 


f d^h d^lL N{Iu...Il) f d^h d^h N{h,...lL) 
J {27r)D ■ ■ ■ {2 t:)D .. 74"'= ^ / (27r)^ ''' (27r)^ 74"i ... 14"'= 



( 2 . 6 ) 


where w is a differential form on V, and contours Cj’s are around the poles of to and also the 
nontrivial cycles of V [29, 31]. Wj are weights of these contours. In particular, to extract 
the coefficients Ci in (2.1), we can find a special set of weights [25-31] such that, 

Ci = '^ ( 2 - 7 ) 

j dCj 


For example, the 474 two-loop massless double box diagram contains 7 propagators. The 
maximal cut is a reducible variety with 6 components [25], each of which is a Riemann 
sphere. The contours are around the intersecting points of these Riemann spheres, which 
are singular points of this variety, and also around the poles of ca. For the 474 two-loop 
double box with massive internal legs, the maximal cut gives one irreducible variety, which 
is an elliptic curve [31]. There is no singular point on this variety, so the contours are 
around the poles of w and also the two fundamental cycles on the elliptic curve. 

Our observation is that if a differential form w on F is integrated to zero, around all 
singular points on F, poles of a; and non-trivial cycles of V. 


f uj = 0, Vj (2.8) 

JCj 

then from (2.7) and (2.4), the integral corresponding to tv can be reduced to integrals 
with fewer propagators. Since from the knowledge of algebraic geometry, it is easy to find 
such w’s satisfying (2.8), we propose a new multi-loop integral reduction method from this 
viewpoint. 
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In this paper, we focus on the cases for which the number of propagators equals DL — 1 
and the maximal unitarity cut gives one irreducible variety. In such a case, the cut solution 
F is a smooth algebraic curve with well defined complex structure. The condition (2.8) 
implies that ui is an exact meromorphic form on V, since the integral 

F{P) = [ uj, yPGV (2.9) 

Jo 

is independent of the path and dP = oj. Then from the study of meromorphic functions on 
V, which is a well-known branch of algebraic geometry, we can list generators for P and 
then derive all forms which satisfy (2.8). 

Explicitly, for this class of diagram, we found that the scalar integral on the cut becomes 
a holomorphic form on V. 


/ 


(27r)^ 


i27r)D D,... D, f 


( 2 . 10 ) 


where the 1-form H is globally holomorphic (without poles) on V. On the cut, the compo¬ 
nents of /j’s become meromorphic functions. We can show that these functions generate all 
meromorphic functions on V. Let P{li, ■ ■ - II) be a polynomial in the components of loop 
momenta, then take the derivative of P, 


dP = fn. (2.11) 

The resulting fQ is an exact meromorphic 1-form. From the analysis above, we get that, 

= 0 -|- (integrals with fewer propagators), (2-12) 

Dk 

so we obtain an integral reduction relation. For the explicit examples in this paper, we can 
show that this method provides all the on-shell part of integral reduction relations. 

In practice, our algorithm can be presented as, 

1. Use integrand reduction method via Grobner basis [40, 41] to rewrite the loop scat¬ 
tering amplitude as the form of integrand basis. The coefficients of integrand basis 
can be determined by fusing tree amplitudes or polynomial division of the Feynman 
integrand. 

2. Calculate the maximal cut of the scalar integral to determine the form of holomorphic 
form n as (2.10). 

3. Calculate the exterior derivatives of all polynomial P’s. In practice, it is sufficient to 
consider the linear P’s and then use the chain rule. Then we get all the on-shell IBP 
relations as (2.11). 

4. Use the obtained integral reduction relations to reduce the integrand basis to master 
integrals. 


f d^h d^h f 

J (27r)^ • • • (27r)^ ... 
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If the integrand contains doubled-propagator integrals, the algorithm will be slightly differ¬ 
ent. We need to solve a polynomial Diophantine equation first, and the procedure will be 
shown in the next section. 

Note that our algorithm is different from the traditional maximal unitarity. Usually, 
maximal unitarity method needs the explicit contour integration to extract the master inte¬ 
gral coefficients. However, our algorithm does not require the explicit contour integration, 
and the explicit form of elliptic/hyperelliptic integral is not needed. The residue compu¬ 
tations to hnd the holomorphic form H and the derivative computations (2.11) are much 
simpler than the contour integrals. 


3 Elliptic Example: Double Box with Internal Masses 

The method explained in the previous section can be used for integral reduction for various 
topologies, for instance, the double-box (Fig. 1) with three different masses for the internal 
propagators. The maximal unitarity of the massless double box was discussed in [25]. Then, 
maximal unitarity structures for double box with 1 ~ 4 massive external legs were studied 
in [26, 28]. In these cases, the unitarity cuts provide reducible curves. 

On the other hand, the unitarity cut of double box with six massive external legs [29] or 
all massive internal legs provides irreducible curves. The integral reduction for symmetric 
double box diagram with internal masses, was discussed in [31], via maximal unitarity and 
the analysis of elliptic functions. Here we show the integral reduction for more generic 
double box diagram with 3 internal mass scales, based on our new method, without using 
elliptic functions/integrals explicitly. 





/C4 


Figure 1. Planar double box diagram with 3 internal mass scales 
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3.1 Maximal unitarity 

The denominators for double box diagrams are 

Di = ll-ml, D 2 = {h - kif - ml, D 3 = {h - h - ^ 2 )^ - mj , 

D 4 = ll - ml, D 5 = {I 2 - kif - ml, Dq = {I 2 - ks- ^ 4 )^ - ml, 

FDj = {h + hf - m 3 . 

We parametrize the loop momenta as, 


- aikl + a2k1 + 

11 =mm 


'2(31) [14] 

and the Jacobian for this parameterization is, 


^2(41) [13]’ 


T 1 9l1 

Ji = det— 

fJ^yi U(Xi 


IS 


4t(s + t) ’ 


_ dll _ 

^ ^ 4f(s + t) ■ 


The solutions for the maximal unitarity cut. 


(3.1) 


(3.2) 


(3.3) 


D\ — D 2 — • • • — Dj — 0. 


(3.4) 


defines an elliptic curve. To see this, we first solve for the variables ai, 02 ) ckS) I3i, 1^2 and 
in terms of 04 and (5^, 


Oi\ — 1, 02 — 0 ) C 13 — 

/?! = 0 ) /?2 = 1) /^s = 

Then the remaining one equation relates 04 and /? 4 , 


mft(s + 1 ) 


045'’ 


m2t(5 + 1 ) 


(3.5) 


iir(o4, (54) = yl(o 4 )/ 3 | + il(o4)/?4 + (7(04) = 0 , 


(3.6) 


Here 71 ( 04 ), ^( 04 ) and ( 7 ( 04 ) are quadratic polynomials of 04 , whose coefficients depend 
on kinematic variables. Formally, /?4 depends on 04 as, 

(«) 

where A is a quartic polynomial in 04 with four distinct roots, for generic kinematics with 
internal masses. Hence the maximal unitarity cut defines an elliptic curve, i.e., algebraic 
curve with genus one, 

C : r/2 = A(o4), (3.8) 

(See the appendix for the basis introduction to elliptic curves.) The choice of keeping 04 
and /?4 and eliminating other variables is purely arbitrary. 
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On the cut, by a short calculation, the scalar double box integral. 


r d% d% 1 

J (27r)4 (27r)4 ... Dj ' 


(3.9) 


has the following structure: 

where the overall factor is not important for the following discussion. 

As [31], it is remarkable that is the only holomorphic one-form associated with the 
elliptic curve. On the cut, the loop-momentum components ai, pi become elliptic functions. 
So we may study the explicit form of these functions, calculate the elliptic integrals and 
perform the integral reduction, in a procedure of [31]. However, in this paper, we propose a 
different procedure: (1) reduce the integrand based on Grobner basis [40, 41] (2) reduce the 
integrand basis to master integrals by the study of differential forms on the elliptic curve. 
The advantage of this approach is that the whole computation is very simple: no explicit 
elliptic parameterization or elliptic integral is needed. The new process can also be easily 
automated on computer algebra systems. 


3.2 Integral reduction 

We now focus on the double box integral with numerator N, 

d% dH2 N 
(27r)4 (27r)4 


m-j 


(3.11) 


Integrand reduction method via Grobner basis method [40, 41] determines that the 
integrand basis contains 32 terms. In terms of (3.2), the basis can be presented as, 

^ otiPl, 04/?!, Pt, alPs, a|, a4/3i a4/3|, alPs, aj, PI,PI, ai/^s, 

04 / 34 , a 4 /? 3 , 04 / 33 , 04 , P 4 ^, P^jO^Ps, a^,a4P4^, a^P^^ 04 , / 34 , P^, 03 , 04 ,1} (3-12) 


On the cut, the integral becomes a meromorphic one-form, 

^'[iVlIcut oc ^^/V(a 3 , 04 ,/ 33 ,/ 34 ) (3.13) 

where is a polynomial in 03 , 04 , P^ and / 34 , and therefore also an elliptic function. We 
now perform the integral reduction, following the strategy as in the previous section. The 
task is to find exact meromorphic one-forms to on this elliptic curves. If two integrals on 
the cut, differ by the contour integrals of such an w, then 

/[Ari]-/[iV2]|,^, = ^a; = 0 (3.14) 

where the second equality holds for all contours, i.e., two fundamental cycles and small 
contours around the poles, because uj is exact. Then the integral reduction between I[Ni] 
and /[A^ 2 ] is achieved at the level of double box diagram. 


I[Ni] — I[N 2 ] = 0 -|- (integrals with < 7 propagators) 


(3.15) 





Note the 04 and /34 generate all elliptic functions on this elliptic curve, as shown in the 
appendix, (A. 5). 

In practice, we find that to find such (u’s, it is sufficient to consider the exterior deriva¬ 
tives of polynomials in 03 , (3^ and / 34 , 


dF{as, a4, 133 , 134 ) 


5F , dF ^ dF OF da4 

a—das + ^— da4 + + ^^«P4 = /- 

oas oa4 ops op4 rj 


(3.16) 


So we need to find the one forms {das,da 4 , d/3s,df34} and then use the chain rule to generate 
integral reduction relations. We can start by calculating da 4 in terms of the holomorphic 
one-form, 

da4 = rj—- = (2A(q!4)/34 -|- B{a4))—- , (3.17) 

T] T] 

where we used the definition rj = \/A and (3.7) to rewrite r/ in function of (34. The purpose 
of this step is to get the a polynomial form of /. 

We can now easily hnd das, 


das 



Qlq 


—Ai— :Tda4 — —-— da4 , Ai = 


mit(s +1) 


a: 


Ai 


(3.18) 


the constant Ai is the product of asa 4 on the cut. To generate the remaining 1-forms, we 
again use the form of elliptic curve. Recall that, 


K(a4, ( 34 ) = A{a4)f3‘l -|- R(a 4)/?4 -|- C{a4) = 0 . (3.19) 


The identity dK = 0 reads, 

d/34 = - (^'(a4)/3| + B'{a4)(34 + C'{a4)) ^ . (3.20) 

Finally we can easily calculate df3s. 




—A2-^d/34 
P 4 


A 2 


(3.21) 


Then use the chain rule, we get all the on-shell IBPs. For example, from (3.16), we analyt¬ 
ically obtain this relation. 


7dbox[al] = ttyTa — 2 -1 ( 3'S^ (mis - mis - - 4mlt + st) Idbox[al] 

2 S^ [4:7712 ~^)\ 

-|- s{4:mls‘^t — 2mls‘^t — 2 m|s^t -|- — 2mlmls‘^ — 2mlmls^ -|- 

—2mlm\s^ + 2mlst^ — Amlst^ — Smlmlst — -|- s^t^)/dbox[«4] 

-|- mlt{s -|- t) (mis — mis — mis — Amlt -|- st) /dbox[l] j + • • • (3.22) 


where ... stands for integrals with fewer than 7 propagators. Consider all polynomials whose 
exterior derivative satisfy the renormalizability conditions, we obtain 23 integral relations. 
Furthermore, consider Levi-Civita insertions which integrate to zero. 


e{l2,ki,k2,ks) h ■ h , e{l2,ki,k2,ks) h ■ k 4 , e(/i, ^2, ^1, ^2), e{li,l2,ki,ks) ■ ( 3 . 23 ) 
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we get 4 more integral relations. So the total number of Mis is, 

#MIdbox = 32 - 23 - 4 = 5 (3.24) 

and explicitly the Mis can be chosen as, 

Midbox — { ddbox [ct4/33] j ^dbox [^^ 4 ] j ^dbox j -^dbox [/^s] j -^db OX [I]}. (3.25) 

or in the conventional choice with X = (li + /2 and Y = {I 2 + /ci)^/ 2 . 

Midbox = {-^dbox[-’i^M],/dbox[-^^],.^dbox[-^],.^dbox[M],/dbox[l]} • (3.26) 

and for instance, the integral reduction in this basis becomes, 

-^dbox[-^^] = 1,^2 -v(^dbox[l] {^m\mls - m\ [ml {s^ + 4st + 16t^) + + 

16s ( 4 m 2 — s) 

+m^s (—m^s + 2m| (m|s + t{s + 4t)) + m\ {2t{s + 2t) — m|s)) + 

—s^t {m\ + m\ (t — 2 m|) + (m| + i))) + 2 /dbox[-^] {mfs (s — 24m2) + 

+2m\ {2m\ (s^ — 2st + 4t^) + s (2m|s — t{s + 2t))) + s (m^s — 2 m 2 (m,|s + 2t{t — s)) + 
+s (m| + Am\t + t^))) — 12/dbox[-^^]s {"ml (s — 8m|) + m^s — 4t) + s (m| + t))) 

+ ... (3.27) 

The whole computation takes about 120 seconds with our Mathematica code. The relations 
are numerically verified by FIRE [3, 4]. 


3.3 Reduction of the double-propagator integrals 


One issue not discussed in [31] is the reduction of integral with internal mass and doubled 
propagators. For the double box diagram, the doubled-propagator integral on the cut also 
becomes meromorphic 1 -forms, so we may carry out the maximal unitarity analysis as that 
in [31]. However, in this section, we show that, our new method is more efficient for reducing 
these integrals. 

Consider the diagram in Fig. 1 with the middle propagator doubled, 


r d% d% 1 

J (27r)4 (27r)4 ’ 


(3.28) 


on the cut, by the degenerate residue computation with transformation law or Bezoutian 
matrix computation [48, 49], we have. 


Idb' 


'OX,217—CRt 


s^w f a4B{a4)da4 


16 


t 


A3/2 


(3.29) 


Unlike the one-forms in the previous subsection, here the one-form have the denominator 
A3/2. It implies that we need to find exact 1-forms like d(F’/A^/^), where F is a polynomial 
in the loop-momenta components. 

Note that A(a 4 ) has four distinct roots, hence A(a 4 ) and A'(a 4 ) have no common root. 
By Bezout’s identity, 


(A(a4), A'(a4)) = (1) 


(3.30) 






and the polynomial Diophantine equation 


/i(a4)A(a4) + /2(a4)A'(a4) = a^B^a^) (3.31) 


has solutions. Such polynomials fi and /2 can be explicitly found by Euclidean division or 
Grobner basis method. The exterior derivative, 


d 



/ 2 A' 

A3/2 


da^ — 2 


JL 

AV2 


da 4 


(3.32) 


determines that, 


Idhox,2\7—cut — 


16 


(/i + 2/2)da4 
AV 2 


(3.33) 


after integrating out the exact form. Now the term A^/^ is removed and we can reduce this 
integral using the result from the previous subsection. In practice, we find a solution such 
that /i + 2/2 is a quadratic polynomial in 04 , so at the level of the double box. 


ddbox,2 — Coddbox[l] T Ci/dbox[-^] T C2ldh OX [A^] + 


(3.34) 


The three coefficients cq, ci and C2 are analytically found by our method and numerically 
verified by FIRE [3, 4]. 


4 Elliptic Example: Sunset Diagram 

The sunset diagram is a two-loop diagram which attracts a lot of research interests [50-72]. 
The sunset diagram with 3 different masses is the simplest loop diagram which cannot be 
expresses in multiple poly logarithms. 

We use our method to study the integral reduction of the sunset diagram (Fig. 2) in 
two dimensional space-time. In this example, we assume that all internal propagators are 
massive. 



Figure 2. Sunset diagram 








Let p be the external momentum, p^ = w?. We can parametrize the loop momenta 
using a variant of Van-Neerven Vermaseren basis [73]. Define two null vectors ei and 62 
such that ef = 0, 62 = 0 and ei • 62 = p^- The Gram matrix of {ei, 62} is, 


In this basis we expand p as, 


and define the auxiliary vector w. 


G= 

\m^ 0 j 

(4.1) 

62 

P=ei + y, 

(4.2) 

tS’]<^ 

1 

II 

3 

(4.3) 


Hence p ■ uj = 0. The two loop momenta can be then generally parametrized as 


— aici + a2e2 , 
h = /3iei + /32e2 , 


(4.4) 


On-shell equations are Di = D2 = D 3 = 0, where Di represent the inverse propagators, 


Di = if - ml, D2 = {p + h- hf - ml, D 3 = . (4.5) 


The on-shell solution can be formally expressed as, 

_ ml _ ml _ -B{a2) ± A{a2) 

^ 2a2P^ ’ ^ 2/32p2 ’ ^ 2 ^( 02 ) 


(4.6) 


where again 02 and (52 satisfy the equation of an elliptic curve A{a2)j5‘l+B{a2)l52+C{a2) = 
0. The discriminant is A = — AAC. 

The sunset integral, on the triple cut, becomes contour integrals of holomorphic 1-forms 

I\3-cut (X. ^ T] = VA(a 2 ) (4.7) 

The integrand basis for the sunset diagram, obtained from Grbbner basis method [13, 
41], contains 12 terms. 


{al,al,ail3i,/3l,a2l5i,a2l32,l5l,ai,a2,l3i,l32, 1} (4.8) 

We now consider the integral reduction for Fig. 2, since the structure is elliptic we 
follow the same strategy as the double box case. As in the double box case, we find the 
one forms {dai, da2, df5i, d/ 32 } and then use the chain rule to generate all the IBP relations. 
We start by calculating da2, 

da2 = -da2 = {2A{a2)(52 + i3(a2))-^ , (4.9) 

r] rj 

we have used (4.6) to rewrite p in function of (32 ■ We can now easily find dai, 

( Ai \ 1 m? 

— I = -Ai—^da2 = -T^da2, Ai = —^ (4-10) 

a2 J Q!2 Ai 2m‘^ 
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Then, 


(4.11) 


dh — — {A!{a2)02 + B'{a2)P2 + C'{a2)) ■ 

Again we can easily calculate djdi 

" (I)=^ 


To generate the 1-forms, given a function F which is a polynomial in ai, 02 ) /^i and (32, we 
can use 


dF = -—doi -— da2 + + ^^«P2 

UCX\ U(y.2 ^Pl ^P2 


(4.13) 


to generate the on-shell part of IBPs. Furthermore, note that the Levi-Civita insertions 


h ■ u, I2 ■ uj , li ■ u h ■ p, I2 ■ u! h ■ p. 


(4.14) 


are integrated to zero. In total, we generate 4 IBPs and 4 independent Levi-Civita insertions 
integral relations. Hence, there are 12 — 4 — 4 = 4 master integrals for the sunset diagram. 
Define that X = li ■ p and Y = I2 ■ P, the four master integrals can be chosen as, 


MI 


sunset 


— {-^sunset [1];-^sunset [^] 7-^sunset ] 7 ^sunset [^] } • 


(4.15) 


For instance, the reduction reads, 

/sunset [^ 1 ^] = ^ (m^ + mlm^ - mlm^ + mlw?) /sunset [ 1 ] 

-I ^ (-3m^ - ml Y ml- mf) /sunset [/^] 

+ 2 ^sunset[/^^] + 2 F mf^ /sunset[l^] + • • • (4.16) 

where ... stands for integrals with fewer than 3 propagators. 

Note that generically, F-dimensional sunset diagrams with 3 distinct masses have 4 
master integral. The four master integrals can be chosen as (4.15) or the scalar integral 
and three doubled-propagator integrals. There is a subtlety that if F = 2, then the 4 master 
integrals are related by Schouten identities [69]. These identities are valid for D < 3, and 
at F = 2 they further reduce the number of master integrals from 4 to 2. 


5 Hyperelliptic Example: Nonplanar Crossed Box with Internal Masses 

We now proceed in studying the integral reduction of the massive nonplanar double box 
(Fig. 3). Unlike the previous examples, this diagram’s maximal unitarity cut provides 
a genus-3 hyperelliptic curve [46, 47]. The structure of holomorphic/meromorphic forms 
on this curve is different from the elliptic case. However, our new approach for integral 
reduction works for this case as well. 

To illustrate our method, we consider the two-loop crossed box diagrams with massless 
external legs and three internel masses scales {mi,m 2 ,m 3 }. Our method also works for 
other crossed box configurations with all massive internal legs. 
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Figure 3. Nonplanar double box 


5.1 Maximal Unitarity and geometric properties 

The denominators for the Fig. 3 are, 

Di=l\-m\, D2 = {h - kif - m\ , D3 = {h - ki - k2f - mf , 

Di = ll-ml, D5 = {I2 - ksf - ml, Dq = {h - I2 + k 4 f - ml, ( 5 . 1 ) 

Dj = (^1 + — ml. 


The on-shell constrains are 


Di 


Dj — 0 , 


(5.2) 


We use the same loop momenta parametrization (3.2). Again, we hrst solve for ai, a2, 03 , 
/?!, 132 and /34 in terms of 04 and /^s, 


Ol\ — 1 , Oi2 — 0 , 

/?! = —(04 -1-03-1-) , 

S 


03 


m\t{s +1) 

3 ’ 

045 '^ 


/52 — 0 , 134 — 


{ml)t{s + t) 


(5.3) 


The rest two variables satisfy a polynomial equation. 


K[a4, (3^) — A{a4)0l “t" B{a.4)f33 + 0 ( 014 ) — 0, (5-4) 


whose solution can be formally represented as. 


/33 


-. 6 ( 04 ) ± A/A(a 4 ) 

2 ^( 04 ) 


A = B^ - AAC 


(5.5) 


Unlike the previous examples, A( 04 ) here is a degree -8 polynomial in 04 with 8 distinct 
roots. Hence the unitarity cut of this diagram provides a genus-3 hyperelliptic curve. (See 
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the appendix for the classification of complex algebraic curves). Note that not all genus-3 
algebraic curves are hyperelliptic, but this one is because of (5.4). 

Before the integral reduction, it is interesting to the see the geometric properties of this 
unitarity cut. Using (5.4) and the statements of appendix, we see that the function 04 is a 
meromorphic function of degree 2 on the curve, i.e., it has two poles Pi, P 2 and two zeros 
Qi, Q 2 - Explicitly, we can check that Qi and Q 2 are distinct, therefore, in the language of 
divisors (A.7), 

( 04 ) = Qi + <52 — Pi — P 2 • (5.6) 


The divisor of 03 is then 

( 03 ) = Pi + P2-Qi-Q2- (5.7) 

The divisor for the function and /34 is more complicated. From (5.4), we determined that 
/?3 on the cut becomes a meromorphic functions of 4 simple poles. The divisor of Ps is: 


m = P2 + Q2 + Wi + W2-Pl-Ql-Zi-Z2. (5.8) 


We hnd that two poles of 04 become a pole and a zero of /? 3 . Similarly, the two zeros of 04 
also become a pole and a zero of The divisor of (^4 is: 


(/34) = Pi + Qi + + ^2 - P 2 - Q 2 - kUi - W 2 . (5.9) 

In summary, there are 8 poles on this hyperelliptic curve from numerators insertions, namely 
Pi, P 2 , Qi, Q 2 , Zi, Z 2 , W\ and Wi. 


5.2 Integral reduction 

First, the integrand reduction via Grobner basis [40, 41] determines that, the integrand 
basis contains 38 terms in the numerator, 

{alPs, al,alf33, a®, a|/33, a|, 04 /?!, a4/3|, a!/?3, al, Pt, Pt, af: “4/31, 

04/31, alh, at /3|, /3|, a4^|, a4^|, at fil, (5.10) 

04/^4, 04/^3, at / 34 , /^s, a 3 , 04 ,1} 


Then, consider the maximal cut for the scalar integral of this diagram. The residue compu¬ 
tation gives, 

Note that unlike the elliptic case, on a genus-3 curve there are three holomorphic 1-forms 
from (A. 13), (which have no pole on the hyperelliptic curve). 


da 4 041/04 a \ da 4 

\/A(a4)’ ■v/A(a4) ’ \/A(a4) 

the scalar integral cut corresponds to the second one, while 


3xbox[a4] 1 7 —cut OC 


041/04 

\/A(a4) ’ 


Ucbox [a3] 1 7—cut OC 


1/04 

\/A(a4) 


(5.12) 


(5.13) 
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It is curious that for a crossed box integral with the numerator linear in li, the maximal 
cut always gives a holomorphic 1 -form. 

This hyperelliptic curve have 6 fundamental cycles and 8 poles as shown in the previous 
subsection. By global residue theorem, only 7 poles’ residues are independent. Therefore 
we may perform maximal unitarity by computing integrals over 6 + 7 = 13 contours. Since 
here we only have one unitarity cut solution, the number of master integers must be less 
than or equal 13. However, in the following discussion, for the purpose of the integral 
reduction, we use our new approach to find exact meromorphic forms instead of calculating 
these integrals explicitly. 

Following what we did for elliptic cases, we would like to generate the IBP relations by 
exact meromorphic 1-forms on the hyperelliptic curve. Again, we calculate the differential 
forms {da 3 , da 4 , d/? 3 , d/ 34 }. 

dai = = (2A(a4)/34 + = {2A{ai)^i - 5 ( 03 ))?^^^^ , (5.14) 

T] rj Ai rj 


where we have used the usual definition rj = \/A- In the second equality, we used the 
on-shell identity, 

Q J- 

(5.15) 


_ mft(s -b t) 
a3Q;4 — M — - n -; 


to recover the form of the scalar integral cut (5.11). The step is not needed for elliptic cases. 
Then, 

das = d ( — ] = —Ai— 5 'da 4 = —^da4 , (5.16) 

V«4/ al Ai 

where again we have used (5.3) to simplify our expression. The exterior derivatives for /3j 
are more complicated. 


d/33 = - (^'(a4)/3| + B'{a4)l3s + C'(a4)) 


and, 


d/34 = d 




as a^da/i 

Ai ri 

m\t{s +1) 


(5.17) 


(5.18) 


Given a polynomial function of {Q;j,/3j}, we can use the chain rule to generate the on-shell 
IBPs. 

We also consider Levi-Civita insertions which are integrated to zero, 


e{l2,k2,ks,k4) l2-h , e{l2,k2,ks,k4) h-ki, e{li,l2,ki,k2), , 

(5.19) 

e{h,l2,ki,ks), e{li,k2,ks,k4), e{l2,k2,ks,k4). 

In total, we generate 25 on-shell IBPs and 6 Levi-Civita insertions identities. Hence there 
are 38 — 25 — 6 = 7 Mis for the non-planar crossed box diagram with three internal mass 
scales. Define that X = {li +p4)‘^/2 and Y = {I 2 +pi)‘^/2, and the Mis can be chosen as: 

{dxbox[-^^], dxbox[i^^], dxbox[ArT], /xbox[-A^], /xbox[-A], /xbox[i^], -^xbox[l]} • (5.20) 
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For this non-planar diagram, the analytic integral reduction relations are significantly more 
complicated. For example, 


^ /xbox (28mfs + 2m\s — 7s^ + — 2m|(s + 2t)) 

J = 8(4m2 - ^ 

3/xbox (4mf (4m| — s — 2t) + s{—2m^ — 2m| + s + 2t)) 

32mf - 8s 

.^xbox (4mf (s + 2t) — s{6m2 — 6rri^ + s + 2t)) 

32mf - 8s ^ 

4xbox [^‘^] (3mf — s + t){28m1s + 277135 — 7s^ + 47713 ^ — 2777 , 2(5 + 2t)) 

16(47773 — 5)5 

~ T 7 ^ a —\ -^-^xbox [ Y ] (47773(35^ + 2st + 4t^) + 5(-27?72S - 277735 - 5t(35 + 2t)+ 

16(4777| — 5)5 

+777|(35^ + 45t + 4t^) + m ^{ — 8 m\s + 35^ + 165t + 4t^)) + 7773(48777|5 + 

+ 5( —35^ + 105t + 4t^) — 27772(5^ + 205t + 8t^) — 2777|(75^ + 205t + 8t^))) + 

+ ——-2-^Txbox [^] (84777,35 + Tlfl\(y —495^ + 245t — 32t^ — 67772(5 + 2t) + 

32(4777| — 5)5 

+ 6777|(5 + 2t)) + 5( —277725 + 5(107773 + 5777|5 + 7 St) + TTr^ —8?77|5 + 55^ + 45t + 8t^)) + 
—m\{s{—7s^ + 345t — 8t^) + 2777|(ll5^ + 8st + 4t^) + 2777^(55^ + 12t^)))+ 


“7T71 —9 -^7xbox[l ](2877735 - 7773 ( 75 ^ + 45t + 32t^ + 27772(5 + 2t) — 2m\{s + 2t))+ 

64(47773 — 5j5 

— 7773 ( 27772 ( 65 ^ + 5t + lOt^) — t{25s‘^ + 485t + 16t^) + 777^(525^ + 465t + 44t^))+ 

— 7773(3277735 + 277725^ + 25t(35^ + 5st + 2t^) — 2?77|(5^ + 165t + 16t^) + TTr^-55^ + 

+385^t + 245t^ + 16t^) + 777|( —135^ + 145^t + 85t^ + 16t^ — 87772 ( 5 ^ + 105t + 4t^))) + 

+ 5(?77|t(10777|5 + 55^ + 2st + 4t^) + 2777|(2777|5 — t(5s + 4t)) + ?772(4?7735 + t(135^ + 

+105t + 4t2) - 27771(35^ + 145t + 4t2)))) + ... (5.21) 

where ... stands for integrals with few than 7 propagators. The integral reduction at the 
level of crossed box takes about 22 minutes with our Mathematica code. We also performed 
the integral reduction of crossed box with doubled propagators, by the same method for 
the double box case as (3.33). All integral reduction relations obtained by our method have 
been numerically verified by FIRE [3, 4]. 

6 Conclusions 

In this paper, we present the relation between the on-shell IBPs and the meromorphic one- 
forms on algebraic curves, for a class of two-loop diagrams: H-dimensional L-loop diagram 
with DL — 1 propagators and one unitarity cut solution. In this case, the unitarity cut has 
a globally well-defined one-dimensional complex structure on it, and hence the analysis of 
IBPs relations can be translated into the analysis of complex curves. 
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By presenting several two-loop examples, planar and non-planar, we show that from 
the knowledge of algebraic curves, it is very easy to construct an IBP relation which reduces 
an arbitrary integral to master integrals. No explicit form of elliptic/hyperelliptic function 
is needed in our method, since only the differential relations of these functions are needed. 
Our method works for the reduction of integrals with or without doubled propagators. 

There are several interesting future directions. In this paper, we mainly consider dia¬ 
grams of DL—1 propagators without three-point massless vertice. If a (DL —l)-propagator 
diagram has three-point massless vertices, generically, the unitarity cut is not an irreducible 
curve but a reducible curve, i.e, union of several irreducible algebraic curves. We find that 
the IBPs obtained from our method, has a smooth massless limit and the limit forms a 
subset of IBPs for these diagrams. For example, the massless limit of our method, applied 
on the massless double box diagram, provides all but 2 IBPs. The missing 2 IBPs contain 
only low-rank numerators, so can be easily found by other algorithms. So in these cases, 
our method would greatly speed up the integral reduction process, even if IBPs are not all 
obtained. In the future, we expect that the algebraic geometry analysis on reducible curves 
will lead the complete set of on-shell IBPs of these diagrams. 

Furthermore, we may consider using geometric properties of algebraic surfaces to study 
loop diagrams with an arbitrary number of propagators. It is well known that the algebraic 
geometry property of surfaces is more complicated than that for curves. However, we expect 
our approach will be generalized for the surfaces cases, because essentially our method does 
not depend on the detailed information of elliptic/hyperelliptic functions or integrals. Only 
the complex structure and differential relations are needed. So the surface cases would be 
studied following this direction, and to recover the “...” terms in our reduction like (3.27) 
and (5.21). 

Finally, we will study the e-dependent part of the integral reduction, based on our 
method. In this paper, we consider diagrams with integer-valued spacetime dimension. The 
ongoing research on two-loop maximal unitarity in dimensional regularization scheme [74], 
also based on algebraic geometry tools, will help us to understand IBPs with dimensional 
regularization from a geometric viewpoint. 
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A Rudiments of Algebraic Curves 

In this appendix we give a brief introduction to the mathematical background of algebraic 
curves used in this paper. The extensive treatment can be found in ref. [33-35]. 

Definition 1 A Riemann surface is a connected one-dimensional complex manifold. 

We are mostly interested in compact Riemann surface. Any compact Riemann surface 
is homeomorphic to a sphere with g > 0 handles attached, and the number g is called 
the genus of the Riemann surface. Since the complex dimension is one, we also denote a 
compact Riemann surface as a complex algebraic curve. However, rigorously speaking, we 
may need to blow up possible singular points on an algebraic curve to make it a Riemann 
surface. 

Definition 2 A holomorphic map between Riemann surfaces X and Y is a continuous map 
f : X ^ Y such that for each holomorphic coordinate (pjj on U containing x on X and ifw 
defined in a neighbourhood of f{x) on Y, the eomposition 

fjwofofPu^ (A.l) 


is holomorphic. 

Definition 3 A meromorphic function f on a Riemann surface X is a holomorphic map 
to the Riemann sphere S' = C U {oo}. 

One very useful theorem regarding to the topological properties of algebraic curves 
is Riemann-Hurwitz theorem. Here we consider the special cases of / : A —)• S. If in a 
neighborhood of P G A, P is located at the origin and / has the expansion f(z) — f{0) ~ z”, 
n > 1, then we say P is a ramified point of / and n is the ramification index of P. Removing 
images of ramification points under /, we get a Riemann sphere excluding a finite number of 
points, namely S. For any point Q £ S, define that d{Q) = #/~^(Q), the number of points 
in the inverse image, d is an integer-valued and continuous function, hence a constant. This 
constant d is called the degree of /. 

Theorem 1 (Riemann-Hurwitz) Let f : X ^ S be a meromorphie funetion of degree d on 
a closed connected Riemann surface A. The ramified points are xi,..., Xn, with multiplicity 
mi,..., m„ . Then the Euler character of X, 

n 

X(A) = 2ci - ^(mfc - 1). (A.2) 

k=l 
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For a compact Riemann surface, the Euler character is related to the genus g, i.e., number 
of handles as, 

X = 2-2g. (A.3) 

The g = 0 compact Riemann surface is a Riemann sphere, while g = 1 compact Rie¬ 
mann surface is an elliptic curve (or torus topologically), g > 1 cases are more complicated, 
and we focus a particular class, hyperelliptic curves, which is defined as an algebraic curve, 

C ■. = h{x) (A.4) 

h is a degree-n polynomial in x with n distinct roots. C has the genus g, if d = 2g + 1 or 
d = 2g + 2, by (A.2). Note that g = 2 curve must be hyperelliptic, but not all 5 > 2 curves 
are hyperelliptic. 

In the hyperelliptic case, x and y become meromorphic function with these properties, 

• X is a meromorphic function of degree 2 on C, 

• y is a meromorphic function of degree n on C, 

• X : C —)• S' has 2y -|- 2 ramified points. If n is even, these points are (x,y) = (ai,0) 
where Oj’s are the roots of h{x). If n is odd, these points are (x,y) = (ai,0) and the 
point at infinity. 

• Every meromorphic function f on C can be uniquely written as 

f = r{x)+ys{x) (A.5) 

where r(x) and s(x) are rational functions of x. 

The last property (A.5) is important for studying the exact meromorphic 1-forms, which 
play the central role of our integral reduction algorithm. 

A.l Riemann-Roch theorem 

We now want to state one of the fundamental theorems of compact Riemann surface X. 
First, we present several definitions. 

Definition 4 A divisor D on a compact Riemann surface X is a formal sum of points with 
multiplicities. 

D=Y, (A- 6 ) 

Pex 

D > 0 if and only if np > 0, VP G X. 

We can naturally associate a divisor to a meromorphic function / in the following way, 

(/)= ^(ordp(/))P. (A.7) 

Pex 


20 - 


where ordp(/) is the leading power of /’s Laurent expansion at P. deg{D) is the degree of 
the divisor defined as 

deg(i:>) = ^ np . (A.8) 

Pex 

C{D) is the space of meromorphic functions / for which, (/) + D > 0. We define 1{D) = 
dim £(!)). Let K be the canonical divisor associated with any meromorphic one form, 

i{D) = l{K - D). (A.9) 

We can now state the theorem: 

Theorem 2 (Riemann-Roch) Let X be a eompact Riemann surface of genus g and D ^ a 
divisor. 

1{D) - i{D) = deg{D) - g + I. (A.IO) 

It is clear that If L> < 0 then 1{D) = 0, so the Riemann-Roch theorem simplihes as, 

i{D) = - deg{D) + g - 1. (A.11) 

On the other hand, if degD > 2g — 2 then i{D) = 0. We have 

1{D) = deg{D) - g + 1 (A.12) 


Theorem 3 If X is a compact Riemann surfaee of genus g then 

1. The space of holomorphic one-forms on X form a finite dimensional veetor spaee of 
complex dimension g, 

2. If ui is a meromorphic differential on a Riemann surface X then the number of zeros 
of u minus the number of poles, counted with multiplieity is 2g — 2. 

For the hyperelliptic curve (A.4), we can hnd an explicit basis of the holomorphic 
one-forms, 

Corollary 1 The g differentials 

x^dx 

-, J = 0,... ,5 - 1, (A.13) 

y 

form a basis of holomorphic differential forms. 

We use this basis frequently in our paper for the integral reduction. 
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